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Abstract. Field-channelled accretion flows occur in a variety of astrophysical objects, including T Tauri stars, magnetic 
cataclysmic variables and X-ray pulsars. We consider a curvilinear coordinate system and derive a general hydrodynamic 
formulation for accretion onto stellar objects confined by a stellar dipole magnetic field. The hydrodynamic equations are 
solved to determine the velocity, density and temperature profiles of the flow. We use accreting magnetic white-dwarf stars 
as an illustrative example of astrophysical applications. Our calculations show that the compressional heating due to the field 
geometry is as important as radiative cooling and gravity in determining the structure of the post-shock flow in accreting white- 
dwarf stars. The generalisation of the formulation to accretion flows channelled by higher-order fields and the applications to 
other astrophysical systems are discussed. 
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1. Introduction 

Accretion is a common phenomenon in astrophysical systems 
ranging from young stellar objects, interacting binaries, galax- 
ies to galaxy clusters. When the magnetic-field stress is larger 
than the ram pressure of the accreting material, the flow is 
confined to follow the magnetic-field lines. The accretion hy- 
drodynamics in these systems are therefore dependent on the 
magnetic-field geometry. 



lar obiects (e.g. iKonigl 1991: Hartmanetal 


]ll994ULill 199(4 


Koldoba et alJl2002HCalvet & Gullbringll998 


[iGullbrineetal. 


2000; iLamzinl 1 19981 Lamzin et all 12001; 


Romanova et al. 


2003; Stelzer & Schmitt 2004), neutron star accretion from 



acting binarie s containing a white dwarf or a neut r on star 
(e.g. lElsner&Lambl Il977t Idhosh & . Lambl Il978l lAronsI 
1993HLovelace et alJl995ULi et alJl996UKrvukov et al.l200fi 
Kold oba et alJ 120021) . In white-dwarf and neutron-star bina- 
ries the magnetic field of the compact star can also be strong 
enough to affect the cooling processes. Despite the fact that 
these stars may have complex magnetic-field structures, the 
dipole field component is important, as it has a longer range 
than the higher-order field components and dominates in the 
regions sufficiently far from the star. 
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Here we investigate the accretion flow onto stellar objects 
in the regime B 2 » %npv 2 (where B is the magnetic field, and p 
and v are the density and velocity of the accreting material re- 
spectively), so that the flow is strictly confined by the magnetic 
field. We apply a curvilinear coordinate system that is natural 
to the dipole-field geometry for the hydrodynamic formation. 
We solve the hydrodynamic equations for an accretion flow in 
which the cooling function has a power-law dependence on the 
temperature and density and thereby obtain the flow-velocity, 
density and temperature profiles. 

We apply the model to the post-shock flow in magnetic 
cataclysmic variables (mCVs), which are close binary systems 
containing a magneti c white dwarf accreting from a red-dwarf 
companion star (see rWarned[l995l) . We calculate the temper- 
ature, velocity and density structure of the post-shock emis- 
sion region and model the X-ray line and continuum emission. 
We co mpare our results to those o btained b y the plane-p arallel 
model Jciievafier & ImamurJl982tlWulll994tlWu et alJll994l 
ICropper et all 19991) which is generally used in spectral analy- 
sis of X-rays emitted from mCVs. 

The paper is organised as follows: In §2 we derive the 
hydrodynamic equations in curvilinear coordinates; in §3 we 
present the treatment of the boundary conditions and discuss 
briefly the numerical scheme that we use; and in §4 to §6 
we show an illustrative example — accretion onto a magnetic 
white dwarf stars — and present hydrodynamic structure and 
spectral calculations. In §7 we consider the use of this formu- 
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Fig. 1. The projection of the (u, w, <p) cordinates on the xz- 
plane. The stellar magnetic dipole is located at the origin and 
the symmetry axis of the dipole is the z axis. The central 
circle respresents the surface of the star. The solid lines are 
the constant-M curves and the dashed lines are the constant- 
w curves. The unit vectors u and w are the unit normals to 
the constant-M and constant-w curves respectively. The <p co- 
ordinate, whose normal is perpendicular to the xz-plane, is not 
shown. 



lation in some other astrophysical applications. Our summary 
and conclusions comprise §8. 

2. Formulation 

2.1. Coordinate systems 

For the analytic study of field-channelled accretion flow, we 
choose a coordinate system with one component along the 
magnetic-field lines. In this representation, the flow is at 
most 2-dimensional (2D); and it can often be reduced to 1- 
dimensional (ID) form when the system has a special geomet- 
ric symmetry. In the present study we consider accretion chan- 
nelled by a dipole field, which is axi-symmetric and, it can be 
shown, essentially ID. We define a curvilinear coordinate sys- 
tem (u, w, <p) for the dipolar field as follows. The first coordinate 
u is defined by the magnetic-field lines generated by a point 
magnetic dipole. The dipole is oriented in the z direction and 
located at the center of the accreting star (the origin of the co- 
ordinate systems), where (x, y, x) are the usual Cartesian coor- 
dinates. As a static magnetic field is the gradient of a potential, 
we use the equipotential curves of the dipole field as second co- 
ordinate w. The third coordinate is the azimuthal coordinate, tp, 
the same as that in the conventional spherical coordinate sys- 
tem (r, 9, <p). The unit vectors u, w and cp are orthogonal to 
each other, as shown in Fig.^ 



In spherical coordinates, the field lines and equipotential 
surfaces of a dipole field satisfy r — ci\ sin 2 6 and r 2 = «2 cos 9 
respectively, where a\ is the radius at which the field line inter- 
cepts the mid-plane of the dipole and «2 is the radius at which 
the equipotential surface meets the polar axis. Thus, we have 

sin20 a m 
u - and (1) 



r 

cos 9 



(2) 



It follows that 9 = cos l (wr 2 ) and the radial location r is a 
function of u and v, as given in (IA.7> in the appendices. Also, 

[ x y]= Vm r 3 [ cos if sin <p J , (3) 
z = wr 3 . (4) 

(See Appendix A for the coordinate transformations and their 
derivations.) 



2.2. Hydrodynamics equations 

The hydrodynamics of the flow is described by the conserva- 
tion equations 



-p+p(V-v) = 0, 
at 

JL + Iv(P e + Pi) = g + f rad 
at p 

dPj Pi d 

-7- -y— -7-p = -(y- i)r ei 

df p at 



d 
df 



P e d 

y— -7-p 

p at 



-(y-i)[v-(q + F rad )-A h -r ei ] 



(5) 



(6) 



(7) 



(8) 



where d/df = d/dt + v • V, and P e , Pi, v and p are the electron 
pressure, ion pressure, velocity and density, y is the adiabatic 
index of the gas, g is the gravitational acceleration, f ra( j is the 
acceleration due to the radiative force, F ra( j is the radiative flux, 
r e i is the electron-ion exchange rate, and Ah is the heating func- 
tion. 

For stationary flows, the time derivatives are zero. For one 
temperature flows, the electron and the ion energy equations 
can be combined by eliminating r e j, yielding 



(v.V)P- — (v.V)p: 
P 



-(7-l)[V-(q+F rad )-Afc] 



(9) 



where P — P e + P{ is the total pressure. We ignore the details of 
radiative transfer for simplicity and assume that the total effect 
of heating and energy loss can be described by a cooling func- 
tion A that depends only on the local hydrodynamical variables. 
We also ignore the radiative force and set f ra( j = 0. Under these 
assumptions, we have a set of three conservation equations 



V • (pv) = , 

(v • V)v + - VP = 

P 

yP 

(v • V)P - —(v • 

P 



V)p = -(r-i)A. 



(10) 

(ii) 

(12) 



This set of hydrodynamic equations is closed after we specify 
an equation of state for the gas. 
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2.3. Dipole-field channelled accretion flow 

As we have assummed by default that the flow follows the field 
lines, the velocity has only the component along the w direc- 
tion, i.e. 



I'W . 



(13) 



If we ignore the effect of stellar rotation, the acceleration is due 
to only the gravitational force, which is radial, i.e. g = -g*r, 
where g* = GM,/R 2 , M* and R„ and the mass and radius of the 
accreting star, and G is the gravitational constant. For accre- 
tion flow along a dipole stellar magnetic field, the gravitational 
acceleration has two components, i.e. 



g = gu U + g w W 



in which 



'* r(u,w) 2 
Qn 



r(u, w) 2 



(14) 



(15) 



(16) 



where r(u, w) is given by equation (IA.7t . The terms Qn and 
<2i2 above are the elements of the transformation matrix Q 
(given in Appendix A). 

The hydrodynamic equations decomposed into the orthog- 
onal components in the (u, w, (^)-coordinate system are 

n 

— (M 3 pv) = 0, (17) 
dw 

v dv 1 BP _ 

hi dw hip dw M 



(18) 



v dP y Pv dp 

r~r ~ i — t = -(r-i)A. 

«2 OW «2 p OW 



(19) 



In the equations above hi, hi and hi are the metric elements 
of the (m, w, (^)-coordinate system, and their explicit expres- 
sion are derived in Appendix A. As we have assumed that the 
flow strictly follows the magnetic-field lines (and omitted an 
explicit consideration of the plasmas and MHD effects, e.g. 
the tension of the field lines and the magnetic stress), we can 
ignore the momentum and force components in u and <p di- 
rections. The three relevant hydrodynamic equations in this 
study are therefore equations dl7> . Jl 8i and ( 1191 . which are 
respectively the mass-continuity, momentum-conservation and 
energy-conservation equations along the magnetic-field lines. 

Direct integration of the mass-continuity equation dl7> 
yields 



hi h^pv — C , 



(20) 



where C is a constant (related to the local accretion rate and the 
geometry of the field) to be determined. It follows that 



d 2 dv 

—(hi h } pv) = hi h-ipv— , 
ow ow 

and the momentum-conservation equation becomes 

d dP 
— (hihipv 2 ) +hih 3 — = hi h 2 h 3 pg w . 
ow ow 



(21) 



(22) 



To s implify the mathematics we consider a variable as in 
ICropper et alJ Jl999l) . which is defined as 



pv 



(23) 



and use it to replace P as a dependent variable in the hydrody- 
namic equations. We also consider a geometrical function 



<H(u,w) = — In (hih 3 ) , 
ow 



(24) 



which describes the variation of the cross-section area of the 
flow-flux tube. (Numerical values of <H for different field 
lines are shown in Appendix B). We may then rewrite the 
momentum-conservation equation as 



dw 



(25) 



Using the mass-continuity and momentum-conservation equa- 
tions in terms of the Iri function and the variable we rewrite 
the energy-conservation equation as 



dv —hi 
dw — v) — v 



r- 1 YH 
——hihiA + - v)v + g w 

C hi 



(26) 



We now have a set of coupled differential equations (equations 
|2~61 and I25> that determine the accretion flow channelled by a 
dipole magnetic field. When the explicit form of the cooling 
function A is given and the boundary conditions are specified, 
equations d25l > and \26\ can be solved along each field line 
(constant u and <p) to obtain the velocity, density and tempera- 
ture profiles. 

2.4. Cooling function 

In this study, we simply consider an ideal-gas law 



pm H 



(27) 



for the equation of state of the gas (where T is the thermal tem- 
perature of the gas, p is the molecular mass of the gas, mg is the 
Hydrogen atomic mass and ks is the Boltzmann constant) and 
ignore the microscopic plasma effects due to the magnetic field. 
In a variety of astrophysical situations, the cooling function of 
an accretion flow can be approximated by 



A = Ap a T p , 



(28) 



where a and are the power-law indices for the density and 
temperature. When a — 2 and [} — 1/2, A is the cooling func- 
tion for the optically thin free-free emission, which is an im- 
portant radiative-loss process in accreting compact stars and in 
some young stellar objects. From the mass-continuity equation 
(I20> and the equation of state for the gas J27> . we can express 
the cooling function in terms of the hydrodynamic variables 
used in our formulation: 



A = A 



C 



h\ /i3 v 



(29) 
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Expressing the cooling function in terms of the hydrodynamic 
variables, the energy conservation is 

(y-lMl 



dv 
dw 



y(£ - v) - v 



h\h-. 



(f 
fH 



vfvP~ a 
v(£ - v) + g K 



(30) 



3. Transonic accretion onto stellar objects 



For non-relativistic, spherical (Bondi-Hoyle type) accretion 
onto stell ar objects, there a re two classes of physical solutions 
(see e.g. Fra nk et al.ll 19921) . In the first class, the flow is sub- 
sonic from infinity to the surface of the accreting star. In the 
second class, the flow is initially subsonic and becomes super- 
sonic at certain radius from the star. As it requires the flow 
velocity to be zero at the stellar surface, the flow will become 
subsonic again via a shock, which converts the kinetic energy 
of the flow to the thermal energy of accreting material. 

The situation is the same for flows channelled by a dipole 
magnetic field, although the dipole field lines, except the two 
originati ng from the poles, a re all closed at finite distances. 
(See e.g. Koldoba et al. 2002J for more details.) In this study 
we consider the second class and focus on determining the den- 
sity and temperature structures of the region between the shock 
and the star. We will show how to construct the shock boundary 
conditions within the hydrodynamic framework that we use and 
solve the hydrodynamic equations with these boundary condi- 
tions. 



3.1. Boundary conditions 

We apply the Rankine-Hugoniot condition to derive the shock 
boundary condition, assuming that the pre-shock flow is super- 
sonic and cold, such that the Mach number M — > oo, we have 

v s = ^y pre , (3i) 

(for an adiabatic index y = 5/3), where V pie and V s are the pre- 
shock and post-shock velocities respectively. As a reasonable 
approximation, we may set V pre = Vg H , where Vff s is the free-fall 
velocity of the gas before encountering the shock. We denote 
the free-fall velocity at the stellar surface as 



V. 



2GM» 



A' 



and at any general position (u, w) Vg 
the boundary condition 

V. 



v s = 



4^ 



(32) 

V*/ yfr. Hence, we have 

(33) 
r s and 



(Here we use the conventions: the shock is located at r 
the stellar surface is located at r = 1.) 

If the collisional energy exchange between electrons and 
ions is efficient compared to radiative cooling then the electrons 
and ions have approximately equal temperatures. The shock 
temperature is given by 



_ 3 fiGM*m H 
s " 8 k B (R* + x s ) 



(34) 



fsee lFranketalJ 1992b where the shock height above the stellar 
surface x s = (r s - l)R t . Recall the definition of the variable ft 
equation J23L and we have 



ft = V s + 



k B T s 



Hm H V s 

which implies that, regardless of the shock location r s , 
ft = Vff s ■ 



(35) 



(36) 



At the stellar surface, we have the stationary-wall boundary 
condition for the velocity, that is v = 0. The variable £ speci- 
fies the specific momentum flux, and determining its bound- 
ary value at the stellar surface is less trivial, as the flow is not 
necessarily perpendicular to the stellar surface. Nevertheless, 
knowing only the value of v at the stellar surface together with 
the values of v and £ at the shock is sufficient to solve the hy- 
drodynamic equations. 

3.2. Modified hydrodynamics equations and numerical 
technique 

The stellar surface corresponds to r — 1 in the spherical co- 
ordinates, but the functional form of the stellar surface in the 
(u, w, <f) coordinates is not trivial. At the surface the indepen- 
dent variable w has different values for different magnetic field 
lines. The boundary value of w for a field line is determined 
(using equation (|5J) only after we specify the colatitude on the 
stellar surface at which the field line is anchored (which we 
shall denote as 9o). 

As the boundary values of w both at the shock and at the 
stellar surface are not easily specified in terms of simple func- 
tions of w and u, we search for an alternative independent vari- 
able. The criterion for the variable is that it increases or de- 
creases montonically in the region of interest. As shown above 
( SD.H the velocity v is better defined at the shock boundary 
and the stellar-surface boundary. It should be smooth in the 
post-shock region and is monotonic in at least in the region 
just beneath the shock and in the region just above the stellar 
surface. 

Without losing qualitative generality, hereafter we will use 
a — 2 and /3 = 1/2 (equation d29> . describing Bremsstrahlung 
cooling) in our calculations. Then the energy-conservation 
equation (I30> can be expressed with v the independent variable 
and w the dependent variable, 



dw 
dv 



= -^[r(£-v)-v]v 3/2 



(y-l)AC 



hih 



L«3 



+ 2^_ v)v 5/2 + ^ v 3/2 



(37) 



Applying the chain rule for differentiation, we can also de- 
rive (from equations 125 \ another equation for the variable ft 
that is 



dv 



dw 
hi — 
ov 

-[r(ft 



v)-v] 



-(f - V) 



^ 2 +^(ft 

«2 



v)v 



3/2 
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(y-DAC 



hih 



|« 3 



fH 

In 



- v) v 5/2 + 



„3/2 



(38) 



The shock-boundary value of v is still not constant among 
the field lines (i.e. for different choices of u). Consider a new 
variable r = v/Vf, where Vf(r) is the local free-fall velocity at 
the radial distance r from the centre of the star, For the shock- 
jump conditions that we assume (eauationsB 1 Iandl33l. the vari- 
able r has a fixed value at the shock and is constant at 1/4. Thus, 
t = at the stellar surface if v = 0. The free-fall velocities 
Vf = V*/ yfr, and hence 



v = 



yHu, w) 



(39) 



where r{u, w) is given in equation (IA.7> . 

The corresponding hydrodynamic equations, for dw/dr and 
d^/dr, can be obtained using the chain rule of differentiation 
with 



dr 

dv 



y/r(u, w) 1 / r' w (u, w)\dw 
tV* 2\ r(u, w) J dv 



(40) 



and r' w (u, w) given in Appendix A. 

The use of r as an independent variable provides a conve- 
nient way to treat the boundary conditions at both the shock 
and the stellar surface. However r (like v) is not guaranteed to 
be monotonic in w for all systems and all choices of u. Thus 
we cannot use it as the independent variable throughout the 
entire post-shock region. In contrast, the coordinate w is al- 
ways monotonic along a field line, but it is not practically us- 
able at the stellar surface, where the velocity gradient dv/dw 
approaches infinity. 

Given that it is straightforward to use w as the independent 
variable in the numerical integration in the entire flow except 
at the stellar surface, we consider a hybrid numerical scheme 
in our calculations. Therefore we use r or v as the variable of 
integration in the vinicity of the stellar surface, < v < dv 
with 6 5 10~ 3 , but we switch to integration in terms of w in the 
upper parts of the post-shock structure. This algorithm allows 
us to solve the flow profile along each field line. 

We calculate the pressure and velocity profiles of the post- 
shock structure by numerically integrating a set of differential 
equations between the upper and lower boundaries. In regions 
where we treat w as the independent variable, we use equations 
d25i and (I26> . Where v or t serves as the independent variable 
then we use equations such as (I37> and J38i or the equivalent 
with derivatives in t. 

Neither boundary condition is completely specified from 
first principles: the position of the shock (expressed in terms 
of w s or r s ) is not initially known. Nor do we have foreknowl- 
edge of the momentum flux at the stellar surface (£,). In prac- 
tice we choose trial values of w s ) and test how the result- 
ing profile conforms to the conditions required at the opposite 
boundary. In one method, we choose trial values of f„ and in- 
tegrate upwards from the stellar surface until a shock is found. 
We adjust until we match the conditions £ s = V,/ yfr\ and 
v s = y[K simultaneously. Alternatively, we may try val- 
ues of r s and integrate downwards from the shock to the point 
where v = 0. If r + 1 when v = then we adjust the trial 



value of r s and retry. These tests are applied iteratively in a 
root-finding routine. 



4. Accretion onto a magnetic white dwarf 

We now apply our formulation to an astrophysical system 
as an illustration. We consider a simple case: accretion onto 
a magnetised white dwarf, in which free-free emission is 
the dominant cooling mechanism. This kind of system can 
be found in intermediate polars and polars, which are cata- 
clysmic variables containing a magnetic white dwarf accret- 
ing material from a low-mass companio n star (for rev iews 
of magnetic c ataclysmic variables, see e.g. Cropper 1990 and 
IWarnerlF9 95i. The post-shock accretion flow in magnetic cat- 
aclysmic variables have been investigated by many workers 
and analytic and semi-analytic results were obtained (e.g jAizu 
1973; Chevalier & Imamura ll982l Iwiilll994t IWu et aljTl994: 



Imamuraetal. 1996; Cropp er et alJl999tlSaxton & WuF l999. 



2001). Their studies address various issues such as the cool- 
ing processes and the effect due to gravity, and the calculations 
either planar or spherical geometries were considered. (For re- 
views of the hydrodynamics of post-shock accretion in mag - 
netic cataclysmic variables, see lWul l2000: Beuermann 2004.) 
Because of the assumed geometry, these studies have not quan- 
tified the effects of the curvature of field lines in determining 
the hydrodynamic structure of shock-heated region. 

For the application of our formulation in the accretion onto 
white dwarfs, we simply set M» = M w( j and R* = R w d, where 
M w d and /? w d and white-dwarf mass and radius, in the hydro- 
dynamic equations and scale the variables accordingly. The 
white-dwarf mass and radius are not independent — when 
the mass is specified one can calculate the radius by means 
of a mass-radius relation. In our calculations, we adopt the 
iNauenbergN 19721) mass-radius relation 



Re 



0.0225 Vl - (M wd /M 3 ) 4 / 3 



(Mwd/Ms) 1 / 3 



(41) 



where M3 = 5.816 M /yU 2 d is the Chandrasekhar mass limit. 
We set the electron mean molecular weight parameter yu w d 
2.00, This is a good approximation for white dwarfs with He, 
C/O or Ne/Mg composition. 

We consider a strong shock for the upper boundary con- 
dition and a cool stationary wall for the lower boundary con- 
dition. The accreting material is an ideal gas with y = 5/ 3 
and approximately solar abundances (Anders & Ebihara Il982t) . 
The mean ionic mass and charge are m = 2365. 8m e and 
Z = 1.0999 respectively, and nJn H = 1.209, Z 1 = 1.3912, 
Z 2 /m = 6.007 x 10 23 g~'. We omit other radiative transport 
and microscopic effect s at the shock and the stellar-surface 
bou ndary (see e.g. Imamur a et al.l Il996t ISaxton & Wul |200 1 ; 
IWu & Cror,nerl200ll). 



4.1. Bremsstrahlung-dominated flows 

If the white dwarf has a weak magnetic field or sufficiently high 
accretion rate then the cooling of the post-shock flow is domi- 
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nated by free-free emission, and Cyclotron cooling is unimpor- 
tant. The cooling function is then 



A = A br = Ap z 



(42) 



with A being a c onstant that depends on th e composition 
of the plasma (see iRvbicki & Lightmanlfl979l For a purely 
Hydrogen plasma, A = 3.97 x 10 16 ,g B in c.g.s. units where 
g B ~ 1 is the Gaunt factor. When we adopt g B = 1.25 and 
approximately solar abundances we have A — 5.06 x 10 16 . 

4.2. Flows with Bremsstrahlung and Cyclotron cooling 

A natural extension of the case for accretion onto weakly mag- 
netic white dw arfs is the ac cretion onto strong-field systems 
(i.e. polars, see Cropper 1990) the flow is cooled by emitting 
both optical Cyclotron radiation and free-free X-rays. We can 
use a composite cooling function to describe the cooling pro- 
cess, i.e. 



A = Ap 2 I — [ 1 + eJ{P,p, u, w) ] 



A to \Pj [p 



h\hi 



(43) 



(44) 



(see lWu et al.ll994l and appendix[Cj, where e s is the parame- 
ter that determines the relative efficiency of Cyclotron cooling, 
and the subscript s denotes quantities evaluated at the shock. 
The penultimate factor in J44l > appears because the Cyclotron 
cooling power depends on the cross-sectional area of the ac- 
cretion stream, which is proportional to h\hi. The last factor 
describes the local magnetic field strength, given by 



B(u, w) 



V4 - 3ur 



(45) 



where B p is its value at the pole. 

Then, we replace p by Clh\hj,v and P/pby — v)v through- 
out the / cooling function, and substitute A into equation ( I26> . 
The resulting hydrodynamic equations can be solved using the 
numerical scheme described in §3. 

4.3. Simplification to previously published models 

We note that the formulation we derive can reduce to the formu- 
lations obtained in the previous studies under certain approx- 
imations and restricted conditions. Take, for instance, equa- 
tions (03 and (OH). If we set hi = h 2 = h 3 = 1, 'H = and g w = 
0, and fix £ to be a constant equal to the free-fall velocity at the 
white-dwarf surface, then equation d38l vanishe s and equation 
( I37l > bec omes t he same as that in lChevalier & Imamural (11982) 
and Wul ill 9941) for the planar flows. If we set u = 0, w = 1/r 2 , 
I12 = -dr/dw, h\hi = I, and 1H — then w e obtain the same 
set of two equations of lCropper et alJ dl999l). 

Comparing the fo r mulati on of ICropper et alJ Jl999l) with 
IChevalier & Imamural Jl982l) and Wu (1994) reveals that the 
former needs two equations to describe the flow while the later 



requires only one. It is because in the presence of gravity the 
quantity which describes the specific momentum flux, is no 
longer conserved along the field lines. In our formulation, de- 
spite the presence of a dipole field, two differential equations 
are sufficient to dete rmine the hydrodynamics, as in the case of 
ICropper et alJ dl999h . We still have two differential equations 
when we set g w = in both equations ( I37t and d38l >. 

To further illustrate the importance of the effects due to the 
field geometry, in the next section we present the numerical re- 
sults of calculations, in which typical parameters of magnetic 
cataclysmic variables are used. We will show that the results 
can differ substantially from those obtained by formulations 
fe.g. lCropperet"aT] [l999) without taking account of the geo- 
metric effects. 



5. Hydrodynamic structure 



The constant C introduced in ( I20t is obtained defining the value 
of pv at the bottom of the accretion column, p*v„. Following 
ICropper et alJ i 19991) we identify the accretion rate m = p*v», 
but in the present work the product pv is not constant in the 
accretion column, since the area of the funnel, oc h\h 3 , varies 
throughout the accretion stream. In our standard illustrative 
cases, we use m = 2.0gcrrT 2 s _1 , a white dwarf mass of 1.0M o , 
escape velocity V* = 6.9 x 10 8 cms _1 and the accretion area is 
taken to be 3.9 x 10 15 cm 2 , exactly 10~ 3 of the global surface 
area. Figure[2]shows how the product pv decreases with height 
within the column for several values of accretion colatitude, 6>o- 

The gravitational acceleration component in the direction 
tangential to the magnetic field lines, which we call g w in (1161 . 
is a fraction of the total gravitational acceleration g. Figure 0] 
shows the relation g w /g. Except for the special case of accretion 
onto the pole, 6q = 0°, the magnetic field lines are nowhere 
vertical to the stellar surface. Thus the component g w is not 
equal to the total acceleration on the surface of the white dwarf. 

Scaled relative to the shock height and parameters of the 
accretion, there are several qualitative differences between the 
post-shock flows in the pres ent dipolar mod el and the cylindri- 
cal accretion model of Croppe r et al.l(ll999l) . 

The distribution of pressure is shown in Fig. |4] The dipo- 
lar accretion model results in proportionally higher pressures 
throughout the post-shock region. 

Corresponding profiles of the density structures are shown 
in Fig. |5] The dipolar funnel results in greater densities 
throughout most of the post-shock region (relative to the 
gas density at the shock), than in the cylindrical accretion 
case. As in previous models with a power-law expression for 
Bremsstrahlung cooling ^e.g. IChevalier & I mamura 1982i lWul 
11994 ICropperetalJ ll999). the density necessarily increases 
asymptotically near the stellar surface. 

Figure |6] shows the gas velocity between the white dwarf 
surface (r = 1) and the shock. The velocity profiles show 
the most significant qualitative difference between dipolar and 
cylindrical accretion models: in the Bremsstrahlung-dominated 
case that we have illustrated, the cylindrical model results in a 
much more constant velocity gradient throughout most of the 
post-shock region. 
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Figure [6] also displays a significant quantitative difference 
between the predicted shock positions in the cylindrical and 
dipolar models. For the cases shown, with accretion colatitudes 
£>o < 18°, the shock occurs at a radius around r s as 1.19. In 
the equivalent cylindrical accretion model, the shock height is 
much lower, r s as 1.105, i.e. just over half the height that we 
would obtain by using the dipole field geometry. 

As the magnetic field strength is increased to some tens 
of MG, the radiative cooling due to Cyclotron emission be- 
comes large compared to the Bremsstrahlung cooling. As e s 
increases, the s hock height decrease s in both the cylindrical ac- 
cretion model (Cropper et alJll999l) and our dipole model. For 
any set of system parameters, (M», m, 9q, e s ), the dipolar model 
predicts the shock to occur higher than in the cylindrical model. 
However the difference is insignificant in high-e s cases where 
the shock height is sufficiently low and curvature effects are 
negligible. For instances in cases with B* = 10, 30, 50MG we 
calculate r s as 1.075, 1.016, 1.007 according to the cylindrical 
model, but the dipolar model yields r s as 1.100, 1.017, 1.007 
when 6 = 0°. 

Figure shows the dependence of shock height upon the 
colatitude of the accretion spot (9q) for several cases of B„ ap- 
propriate for polars. The maximum shock height generally oc- 
curs for accretion along the field line with #o = 0°. In cases with 
low B t the variation of r s with 0q is less than in the stronger field 
cases. When the magnetic field is stronger the shock height is 
smaller compared to the white dwarf radius, and the dipolar 
funnel effects less significant. 

Figure [8] shows equivalent relations between shock height 
and 6*o, but for a greater accretion rate, m = 5.0 gem -2 s _1 . 
Bremsstrahlung cooling is more efficient in this denser flow, 
and the shock heights a re lower than in Fig.|7] as expected from 
previous studies (e.g. lAizulll973l IWuHl994l Emamura et alJ 
Il996l) . Again, as in Fig.0 r s is more variable in 6q in the cases 
with weaker magnetic field and greater shock height. 

Figure [9] shows the temperature as a function of height 
within the column in our standard illustrative case. For compar- 
ison we calculate the temperature profile with the same param- 
eters for the white dwarf and accretion rate, but applied in the 
cylindrical accretion model. In this specific case the shock tem- 
perature in the cylindrical accretion model is 4 keV greater than 
it is in the dipolar accretion model. The dipolar accretion model 
generally yields a lower shock temperature than cylindrical ac- 
cretion, because of the greater shock height and thus lower free- 
fall velocity at the shock. The temperature predictions of cylin- 
drical and dipolar accretion models are much more alike for 
cases where the shock is closer to the white dwarf surface (and 
dipolar funnel effects are less significant). 

Figure EH compares the cylindrical and dipolar models' 
results for the temperature structure of the accretion column, 
for four different values of the white dwarf mass, in the 
Bremsstrahlung-dominated limit. The dipolar model generally 
has a lower shock temperature than the equivalent cylindrical 
model. The difference between models is greatest for the most 
massive white dwarfs, as their accretion shocks are higher (and 
thus more affected by funnel geometry). 

In the next section we present calculations of the X-ray 
spectrum emitted from the post-shock flow, for two represe- 



M t =^.0 Mq , m = 2.0 g cm~ 2 s~ 1 




0.0 0.2 0.4 0.6 0.8 1.0 

(r-1) / (r s -1) 



Fig. 2. Variation of the local mass accretion rate pv within the 
post-shock structure. Different cases of the accretion colatitude, 
6 Q = 0°, 18°, 36°, are marked by full, dashed and dotted lines 
respectively. At the white dwarf surface we have p»v* = m by 
definition. In this illustrative case m = 2.0 gcrrT 2 s _1 , B, - 0. 
and M» = 1.0M Q . The horizontal axis is the altitude above the 
white dwarf surface, relative to that of the shock, r s - 1 . 



native cases of M*. We compare the predictions of the dipolar 
accretion model with those of the planar models, thereby show- 
ing the precise effects of curvature. 

5.1. Comments on curvature effects upon the flow 

In section §[5]we use an accreting magnetic white-dwarf star as 
an illustration and demonstrate that the role of the field geome- 
try can be as important as (or even more than) the role of gravity 
in determining the velocity, density and temperature profiles of 
the flow. We can elaborate this more clearly by inspecting the 
meanings of the terms in the hydrdodynamic equations. Each 
of the terms in the square brackets in the right side of equa- 
tion J26l > corresponds to an energy-transport process. The first 
represents the radiative loss, the second is determined by the 
geometry of flow, and the third is due to an external force field, 
which is only gravity here. A characteristic of the second term 
is the H function. In the curvilinear coordinates that we use, hi 
and /z3 are the metric elements corresponding to the coordinates 
perpendicular to w. Thus 1ri{- d\s\(h\h{) j 'dw) is the measure 
of the changes in cross-section area of the magnetic-flux tube 
along the flow, and it plays an important role in determining the 
efficiency of compressional heating. 

It is worth noting that the 'H function also equals to h^S • b, 
where b is the unit vector tangential to the magnetic field, and 
I12 is the metric for the coordinate component along the field 
line. We may also use this relation to determine the H function. 
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Fig. 3. Variation of the component of the gravitational acceler- 
ation tangential to the magnetic field line. The five curves corre- 
spond to cases of hot-spot colatitude 8 Q = 0°,6°, 12°, 18°, 24°. 
The parameters (M„ B*, m) take the values used in Fig.|2 



For accretion onto stellar objects, g w ~ and £ ~ Vs, 
where Vff is the mean free-fall velocity. The metric element 
h\ ~ 0(1), and hence, the geometry function fi ~ 0(1). From 
equation J38I . we can deduce that the effects due to field cur- 
vature (geometry) and gravity are often comparable. Generally, 
the geometry effect is more important for regions close to the 
stellar surface than for regions further away, because the con- 
vergence of the dipole field lines. 

We note that equations J25l > and J26l > are general equations, 
in principle, applicable for accretion flows channelled by any 
field geometries. For example, if the magnetic field has pla- 
nar parallel structure, then we have hi = /12 = hj = 1 . It fol- 
lows that 'H — 0. In this flow, the accreting material will not 
be compressed by the magnetic field. If the flow is spherical, 
then hi = 1, h% = r, hi — rsind (see e.g. lArfken & W eeer 
and H = 2/r. The flows channelled by dipole field have 
been shown in the sections above. The generalisation to chan- 
nelled flows in higher-order fields is possible. What one needs 
are first to find the functional expressions of the field lines and 
the equipotential surfaces, and use them to define the coordi- 
nates. After the coordinate system is specified, we can derive 
the metric elements h\, I12 and /13 and determine the <H func- 
tion. 

The effects of more complex field geometries depend on 
the alignment and distribution of the higher-order field com- 
ponents. If the field lines near the accretion spot diverge more 
rapidly than in the dipolar model (i.e. a more flared accretion 
funnel) then we expect an accentuation of the effects found in 
the present work: a higher shock position and lower shock tem- 
perature. If the field lines near the accretion spot diverge less 



Fig. 4. Post-shock pressure profiles with the parameters 
(M,, B„, m) used in Fig. [5] The vertical axis is scaled in terms of 
the free-fall velocity (Vs s ) and density (p s ). immediately down- 
stream of the shock. The lower curve (CWRK) is the cylindri- 
cal accretion model. The upper curve is the dipole model with 
80 = 0°. Varying 6q has little effect on the appearance of this 
curve. 

M«=1.0 Mq , m = 2.0 g crrr 2 s _1 
1 00 : 'J 




0.0001 0.0010 0.0100 0.1000 1.0000 
(r-1) / (r.-l) 

Fig. 5. Density profiles of the accretion columns shown in 
Fig. 0] scaled to the gas density immediately downstream of 
the shock (p s ). In the upper curve 6q = 0°, and variation of 80 
has little effect. The cylindrical accretion model (CWRK) is the 
lower (dot-dot-dot-dash) curve. 
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1.20 



Fig. 6. Velocity profile of the post-shock structure, normalised 
to the escape velocity and the stellar surface, V,. The upper 
curve represents the model of Cropper et al. (1999). The lower 
curve shows the effects of accretion in a dipole-field accretion 
funnel, with 9q = 0°. The system parameters are the same as in 
Fig.E 



rapidly then the accretion flow will be a closer approximation 
to the older, planar accretion models. 
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6. X-ray emission 

6.1. Spectral calculations 

We use the analytic model to provide electron density and tem- 
perature values throughout the post-shock flow. The structure is 
divided into between 2 x 10 3 and 2 x 10 4 strata, depending on 
the adaptive integration steps chosen by the routine that solves 
the set of differential equations. Using these (« e , T e ) values, an 
XSPEC subroutine for a MEKA L model of an optically thin 
plasma jMewe et alJll995t IPhillips et al.lll999l) calculates the 
X-ray emission spectrum of each layer. The total X-ray spec- 
trum of the accretion column is a sum over all the layers' spec- 
tra, weighted by cell volume. We omit the few cells for which 
n e > 1.0 x 10 18 cm" 3 , for which the MEKAL model is inappli- 
cable and fails. In reality the gas at these densities is optically 
thick and merges into the stellar atmosphere. 

Figure ^2 presents the X-ray spectrum calculated for a 
high-mass white dwarf (l.OM ) with a magnetic field strength 
of 10 MG at the accretion hot-spot. The uppermost panel rep- 
resents the dipolar accretion model (I). The lower panels are 
ratio plots, comparin g the dipolar model with: (II) the model of 
ICropper et alJJl999 ) (planar geometry, including gravity); (III) 
the model oflWuNl 994 ) (planar, but without gravity); and (IV) 
an isothermal, homogeneous post-shock model. For photon en- 
ergies ranging from 0.2 to 10 keV, the spectrum of model (II) is 



Fig. 7. Effect of the accretion colatitude, Go, upon the ra- 
dial location of the shock, r s (in units of stellar radius), for a 
white dwarf of mass M„ = l.OM and m = 2.0 gem -2 s _1 . 
Results are shown for cases with magnetic field strength B„ = 
0, 10, 30, 50 MG at the accretion spot. 

harder than that of the dipolar model (I), and has less line emis- 
sion. It may be unsurprising that the CWRK model has a harder 
spectrum, since the shock temperature is higher than in the 
dipolar accretion model. The planar accretion model without 
gravity (III) produces a softer spectrum than the dipolar accre- 
tion model, and more line emission. The homogeneous shock 
model (IV) yields a spectrum that is qualitatively very different 
from any of the inhomogeneous shock models. It shows a much 
harder spectal slope. 

Curvature effects are strongest for those cases where the 
shock height is significant compared to the stellar radius, and 
thus the differences between spectra in comparable CWRK and 
dipolar accretion models is most noticeable for the cases of the 
most massive white dwarfs. For lower mass white dwarfs the 
shock height is small compared to the stellar radius, curvature 
effects are minor, and the dipolar and CWRK models predict 
similar X-ray spectra. For the case of M* = O.5M (Fig. I12> 
the spectra predicted from these models are almost indistin- 
guishable. 
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M«=1.0 Me , m = 2.0 g cm -2 s" 



Fig. 8. Radial location of the shock as a function of 6>o, with the 
same conditions as in Fig.0 except the accretion rate is higher: 
m = 5.0 gcnT 



,-2 s -i 



6.2. Effects on X-ray spectrum 

For white dwarfs of low or medium mass (e.g. 0.5M o ) the 
shock temperatures, post-shock structures and X-ray spectra 
are nearly the same in the dip olar accretion model a s in the ear- 
lier planar accretion model of Croppe ret alJ Jl999h . For white 
dwarfs of higher mass (e.g. 1 .OM ), the dipolar accretion model 
predicts greater shock heights, lower shock temperatures and 
softer X-ray spectra in the band from 0.2 - 1.0 keV. This im- 
plies that for a given observed spectrum, the mass estimate M* 
is necessarily greater when dipole field curvature is taken into 
account. 

The lower shock temperature does not at first sight provide 
a route to resolving the high white dwarf masses derived using 
spectral fits to X -ray data, and which are the subject of some 
controversy (e.g.fSchwope et al. 2002), since even higher mass 
white dwarfs are required to provide sufficient hard X-ray flux 
to fit the observed spectrum. 



60 



50 r 



CWRK 




0.4 0.6 
(r-1) / (r s -1) 



1.0 



Fig. 9. Variation of the gas temperature in the accretion column 
as a function of altitude scaled to the shock location. System 
parameters (M t ,B*,ih) are the same as in Fig. [2] The solid 
curve is the case of 6q = 0°; results for 8q < 18° appear in- 
distinguishable on this plot. The nearby dotted curve is the ex- 
treme case of 6 — 60°. The upper curve (dot-dot-dot-dashed) 
portrays results from the cylindrical accretion model (Cropper 
et al. 1999). 
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Fig. 10. Comparison between cylindrical and dipolar geome- 
try for four different values of white dwarf mass. From bot- 
tom to top each pair of curves represents: M* = 0.7,0.8,0.9 
and 1 .QM Q . In each pair the upper (dotted) curve represents the 
cylindrical model (CWRK) and the lower (full line) curve rep- 
resents the dipolar model. We take 9q = 18°. 
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Fig. 11. Comparison of X-ray spectra calculated using (I) 
the dipole accretion model with gravity; (II) CWRK accre- 
tion model with gravity but without curvature effects; (III) pla- 
nar accretion model of Wu (1994) without gravity or curva- 
ture; (IV) isothermal model using the post-shock temperature 
and density derived from model (I). The upper panel shows 
the X-ray spectrum of model (I) directly. The remaining pan- 
els present ratios between the dipole model spectrum and the 
other models. In all cases we set m = 2.0 g cirT 2 s~ l ,B t = 
10MG,M« = 1.0M o ,6> = 0°. 

7. Application to Other Systems 

The formulation that we derived is not restricted to accretion 
in magnetic white dwarfs. It is applicable to a variety of astro- 
physical systems, provided that the flow is strictly confined by 
the magnetic field lines, and the radiative heating and cooling 
processes can be parametrised in terms of the hydrodynamic 
variables. 

An example is accretion in young stellar objects. There are 
strong observational evidence that accretion flow in T Tauri 
stars are chann elled by the stellar magnetic fields (e.g. 
Basri etal Jl992h . The main differences between the mCV case 
and the accretion in T Tau stars are: that the cooling process 



Fig. 12. Comparison of X-ray spectra calculated according to 
the four models, labelled as in Fig. ^2 but with stellar mass 
M, = O.5M . 



is dominated by line cooling instead of free-free and cyclotron 
cooling. Moreover, irradiative heating of both pre- and post- 
shock flow is important in determining the flow hydrodynam- 
ics (see e.g. Martin 1996). Under these conditions, we need to 
relax the assumption of a completely dissociated gas and gen- 
eralise the hydrodynamical equations to multi-ion and electron 
flow. In a fully self-consistent treatment the equations of hydro- 
dynamics need to be solved simultaneously with the ionization 
structure equations (e.g. lLamzinir9 98). 

Another example of application is the accretion onto a 
slowly rotating neutron star in a binary system (i.e. X-ray pul- 
sars). The accretion flow is magnetically funnelled, from ei- 
ther the wind of the c o mpanion star or an acc r etion disk (e.g. 
Eisn er & Lambl 1 19771 iGhosh & Lambl Il978l Lovelace et all 
1995; Koldoba et al.l l2002l) . In cases of disk accretion, the gas 
is funnelled from the inner edge of the disk, and the radius 
of this location effectively determines the colatitude of the ac- 
cretion hot-spot on the stellar surface (9q in our formulation). 
For details of interaction of the channelled flow and the accre- 
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tion disk in X-ray pu l sars, see e.g. | Anzer & BoerneJ £983); 
ISoruit & Taaml ( ll990l) : lliet alJJl996l) . 

8. Summary 

We have investigated accretion onto stellar systems in which 
the flow is strictly channelled by a dipole magnetic field. Such 
flow occurs in regions close to the stellar surface, where B » 
8npv 2 . We derive a set of hydrodynamic equations using the 
curvilinear coordinates natural to the field geometry. The equa- 
tions are solved to determine the flow velocity, density and tem- 
perature profiles. We show that the dipole-field geometry can 
cause significant compressional heating, and this effect can be 
comparable to radiative cooling and gravity in determining the 
structures of the flow near the surface of the accreting star. 

The formulation that we have derived is applicable to a va- 
riety of astrophysical systems, from white-dwarf stars in mag- 
netic cataclysmic variables to young stellar objects. It would 
be possible to generalise the formulation for dipole-field chan- 
nelled flow to the flow channelled by a higher-order field using 
curvilinear coordinates. 

We demonstrate that our analytic model can be efficiently 
interfaced with spectral software to reproduce fits in analysis 
of high quality X-ray spectra, obtained by satellites such as 
Chandra and XMM-Newton. This is a next-generation model 
for fitting accretion parameters, including white dwarf mass, 
magnetic field strength, accretion colatitude and mass flux, to 
specific observed systems. 
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Appendix A: Orthogonal Basis and Metric 
Elements 

Let i, j and k be the standard orthogonal unit vectors in the 
Cartesian coordinates. The unit vectors of our dipolar curvilin- 
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ear coordinate system (u, w, ip) and {%, j, k) are related to the 
Cartesian unit vectors via 



(A.l) 



where (see e.g. lArfken & WereihoOlh the transformation ma- 
trix is 



u 




i 


w 


= U 


3 






k 



The first derivatives of the radius function, r' u = dr/du and 
r' w = dr/dw, are obtainable by implicit differentiation of iA.6i : 



and 



4w 2 r 3 + it 
-2wr 4 



4w 2 r 3 



(A. 12) 
(A. 13) 



U = 



±dx J_dy _\_dz 

h\ du h\ du h\ du 

±dx ±dy_ _\_dz 

h.2 dw hi dw hi dw 

J_dx 1 dy J_dz 

hj ~Bip h?, ~d~ip hi ~Bip 



The relation between (x, y, z) and (u, w, p) is given by equa- 
tions Q and @. It follows that 



<A " 2) h[-y\-\&+*>[ 



cos ip sin ip 



The coefficients h\, hi and hj, are the metric of the curvilinear 
coordinates system, and they are 



dz 
du 



i 2 i 

= 3wrr' 



dx 



dy 



dz 



,(A.3) 



cos ip sin ip 



d[u,w,p]J \d[u,w, <p] J \d[u, w, ip]/ p. 

-£ = r\r + 3wr' w ), 
The transformation from between dipolar coordinates and dw 



spherical coordinates (r, 0, <p) is given by: 



r 




i 




u 





= N 


3 


= Q 


w 


<P 




k 








— [ x y z ] = Vwr 3 [ - sin ip cos <p J 



(A. 14) 
(A. 15) 
(A.16) 
(A. 17) 
(A.18) 



(A.4) 



where the preceding matrix is 
N 



Consider these substitions into (IA.3t : (i) equations (IA.14> 
and (IA.15> into h\ ; (ii) equation 1A. 161 and ( IA. 17l > into h%\ and 
(iii) equation jA.18> for hy, Then we obtain metric elements 
explicitly: 



sin ft cos tp sin ft sin p cos 6 
cos ft cos p cos ft sin <p - sin 
- sin ip cos ip 



(A.5) hx = JJL( r + 3 u r' u ) 2 + (3wr 2 r' u ) 2 , 



Q = N • M, and M = U' = IT 1 (the transpose of U), Thus Q ^ _ C 7 r + 3 w r > )2 + l urr '2^ 
completely prescribes the relationship between (u, w, ip) and 1 4 

With m and w defined according to Q and we have 
sin 2 ft = ur and cos 2 ft = w 2 r 4 , and as cos 2 ft + sin 2 6—1, 
we therefore have 



— \u r 3 . 



(A.19) 

(A.20) 
(A.21) 



It is also worth noting that it, w and cp vary along the cor- 
responding coordinate axes. It can be shown that 



w 2 r 4 + u r - 1 = . 



This is a quartic equation, which has two real and symmetri- 
cal ro ots and two imaginary roots (e.g. lAbramowitz & St egun 
11964 . The real root that is of physical relevance is 



(A.6) 5fl 
du 



w dh\ ip dhA 
I12 dw /13 dp j ' 



r= y/Y-S + 
where 



4w 2 Y 



Y = y/(S + +S.) 2 + 3(S + -S-) 2 



(A.7) 

(A.8) 
(A.9) 



dw 

dw 

d<p 
dip 

(see e.g. lArfken & Wegerl200l1) . 

Appendix B: The Function 



u dli2 (p dll2 

h\ du hi dip 

u dh-$ w dhA 

hi du h2 dw j 



(A.22) 
(A.23) 
(A.24) 



X 



64w 4 



128w 4 



and 



M 4 + 



The function 1-l(u, w) is a purely geometrical function, depen- 
(A 10) d ent on tne curvature of the coordinate system. In the (m, w, tp) 
coordinate systems that we use, it is 



256w 2 
27 



(A.ll) <H = 



3ur 5 
(hih 3 ) 2 



(B.l) 
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Appendix C: Cyclotron cooling efficiency £ s 




r radial coordinate 



Fig. B.l. From top to bottom, the < H(u,w) function ( for 
6 = 0°,6°, 12°, 18°, 24°, 30° and 36° respectively. The hori- 
zontal axis is the radial distance from the centre of the star. In 
the calculation m - 2.0gcirT 2 s~ 1 andM* = 1.0 M Q . 



where 



1 2 r' r" 
f — _|_ l j_ _|_ liw 

u r r'„ ' 



1 7 r' r" 
— __ -l_ _ _ _)_ uw 

w 2 r r' 



(B.2) 



(B.3) 



In the above expression, r" w = dr'Jdw = d(dr/du)/dw, which 
is obtained by implicit differentiation of jA.12l i: 



[<, + (8wr 3 + 12wV<,)r; ( ] 



4w 2 r 3 + // 



(B.4) 



We can also relate the change in p v to the 1-i(u, w) function. 
The mass-continuity equation implies 



d d 
hihj, — {pv) +pv— (hih 3 ) = 
ow ow 



Hence, 
d 

—ln(pv) 
Ow 



-<H . 



(B.5) 



(B.6) 



Figure i lB. 1 1 shows the "H(u, w) function for some values of 6„, 
the angle between the axis of the dipole and the position vector 
of the foot point of the magnetic-field line. 



In our illustrative system the Cyclotron emission from the post- 
shock region is optically thick at frequencies up to some cut- 
off, a>, (e.g. IChanmugam & Wapierll 19791 IWadaet alj fl 980: 
lLanger et alJl982l). As derive d in lsaxtonl i ll 9991) and theappen- 
dices o nCropper et alJ Jl999l) . the effective volumetric cooling 
function is 



A, 



ttD k B T 

a I2n 2 c 2 



(CI) 



where cj, * 9.87w c (©/10') uui (r/10 s K) us , cj c = eB/m e c is 
the Cyclotron frequency, ® = 27ren e D/B, D is the diameter of 
the accretion flow, a = 7rD 2 /4 is the cross-sectional area, B 
is the magnetic field strength, T and « e are the electron num- 
ber density and temperature respectively. Calculating A cy in the 
conditions immediately downstream of the shock, and taking 
its ratio with the Bremsstrahlung cooling ( |42> gives the relative 
efficiency parameter e s = A cy s /Ab ri . s 1994). An appropri- 
ate series of substitutions leads to ( 1441 . 

However there is a practical problem to performing cal- 
culations starting from a chosen value of e s , since J44> de- 
pends on functions of the shock position, which is initiall y un- 
known. In the earliest cylindrical accretion models jAizull973t 
IChevalier & Imamurall 983 IWvJ 19941) it was possible to elim- 
inate the shock height by normalisation, but this is impossi- 
ble for formulations that include explicitly position-dependent 
effects such as gravity JCropperetalJ fl999). For computa- 
tional convenience, we modify the definition by replacing the s- 
subscripted quantities in J44l > with reference values at the stellar 
surface /i3»), and ideal post-shock pressur e and densi - 

ties (P' s ,p' s ) of the zero-gravity, cylindrical model o f tWul(ll994 . 
In that model the shock was assumed to occur close to the stel- 
lar surface. The corresponding efficiency parameter is denoted 
e s „, and it yields, after solving for the actual shock location, the 
value of 



4 - 3wr s 
(4 - 3w>f 



(C2) 



In our main illustrative case, with M* = l.OM , V t = 
6.910 x 10 s cms -1 , m = 2.0 gcm^s" 1 and approximately 
solar abundances, we derive e s * = 3.2,72.7,312 and e s 
1 .6, 64.6, 297 for B« = 10, 30, 50 MG respectively. 



